Potential flows in the Reissner-Nordstr6m-(anti) de Sitter metric: numerical results 
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Numerical solutions for the integral curves of the velocity field (streamlines), the density contours, 
and the accretion rate of a steady-state flow of an ideal fluid with p = Kn~' equation of state are 
presented. The streamlines and velocity flelds associated with a black hole and a rigid sphere in a 
Reissner-Nordstr6m-(anti) de Sitter spacetimes are studied in some detail. For each case the fluid 
density is studied using contour lines. For 7 7^ 2, we found that the studied properties of the 
fluid are more sensitive to variations of the electric charge and the cosmological constant. Also, 
the accretion rate was found to increase or decrease when the cosmological constant increases or 
decreases respectively. 
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I. INTRODUCTION 

Most of the articles about potential flows in this area 
deal with the important case of fluid motion evolving in 
the spacetime associated with compact stars and black 
holes. The implementation of new background metrics 
brings some new challenges. First, metrics other than 
Schwarzschild and Kerr are not so well studied, some- 
times a complete understanding of the physical meaning 
of the metric is missing. Also, the solutions of the fluid 
equations in a non trivial metric may be quite involved. 
In particular, the search for significant boundary values 
(or initial conditions) presents a non trivial problem. A 
simple and paradigmatic case of a flow is the stationary- 
zero-vorticity flow of a fluid with adiabatic stiff equa- 
tion of state. In this case, for relativistic flows, the fluid 
equations admit analytical solutions for some particular 
metrics [H 0, Hi- These solutions are used as test-beds 
for testing almost all the numerical hydrodynamic codes 
in the subject. Other potential flows in a nonstation- 
ary background and different equation of state had been 
studied, see for instance [i, £]■ Also, the solutions for 
potential flows permit to test new optimized codes in re- 
solving nonlinear hyperbolic systems of conservation laws 
^, see for a representative sample of numerical schemes 

In previous articles [E 0, H, Q the authors used for 
the fluid the stiff equation of state, i.e., a polytropic, 
p = Kri* with 7 = 2. This stiff equation of state over 
simplifies the partial differential equation for the veloc- 
ity potential field $ (it becomes linear). In this case 
the sound velocity in the fluid is equal to the speed of 
light. Therefore the stiff equation of state represents a 
limit situation not easily encountered in usual physical 
situations. Some articles [1, Q studied the more realistic 
situation, a fluid with a polytropic equation of state with 
1 < 7 < 2 in the presence of either a rigid sphere or a 
black hole. The sound velocity in this case is less than 



the speed of light and the partial differential equation for 
$ turns to be nonlinear and needs to be solved numer- 
ically. In this work we present a nontrivial extension of 
the investigations about potential flows by studying the 
streamlines and baryon density contours for a fluid with 
polytropic equation of state with 1 < 7 < 2 in the pres- 
ence of either a rigid sphere or a black hole in a Reissner- 
Nordstr6m-(anti) de Sitter background. We consider also 
the accretion rate of particles into the black hole and its 
dependence with 7, the electric charge and the cosmolog- 
ical constant. We assume that the fluid is a test fluid, i.e., 
the metric does not evolve and it is given a priori. The 
state equation and the idea of a rigid star (rigid sphere) 
we use are idealized. However, they bring important re- 
sults about the behavior of the fluid, and the difficulties 
involved in this kind of scenarios. 

This work is divided as follows. In Sec. II we present 
the basic equations that describe potential flows and the 
nonlinear equation for the velocity potential for a poly- 
tropic equation of state of the form p — Kri^ . In Sec. Ill 
we present the numerical method used to solve the non- 
linear partial differential equation presented in Sec. II. 
In Sec. IV we show the numerical results for the density 
contours, the streamlines and the accretion rate for the 
black hole and rigid sphere cases for different values of 
the parameters A, M, Q and 7. Finally, in Sec. V we 
summarize our results. 



II. BASIC EQUATIONS 



We start from the energy-momentum tensor for an 
ideal fluid, T^^ = (p + pW^JJ^ + 

P9tJ,i^j where p is the 
pressure, represents the fluid four- velocity, p — po + e 
is the total energy density, pQ is the rest mass energy 
density, and e is the internal energy density. Our conven- 
tions are G — c — 1, metric with signature +2. Partial 
and covariant derivatives are denoted by commas and 
semicolons, respectively. The conservation equations, 
T.'^" = 0, for this kind of fluid reduce to 
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ip + p)U^, + p,,,U^ = 0, 



(1) 



2 



and 



(2) 



which are, respectively, the 'mass' conservation law 
and the Euler equation. For isentropic flows we have 
(cr/n);^ = 0, where a is the entropy per unit volume and 
n is the baryon number density. In this case the equa- 
tions of motion ([2]) take the form Q 



U^'uj., = 0, 



(3) 
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where w^i, is the relativistic vorticity tensor defined as 



(4) 

The potential flow solution of this equation {ui^i, — 0) is 

hU^ = (5) 

where h — [p + p)/n is the enthalpy per baryon. From 
Eq. ([5]) and the baryon number density conservation 
equation, 



(nU^),^ = 0, 



(6) 



we obtain the differential equation for the velocity poten- 
tial, 



□$+[ln(^)],„$'"-0. 



(7) 



where, □$ = [^y^g'^''^^^]^v/^y^^. This equation is 
in general nonlinear and depends on the fluid equation 
of state. If we consider a polytropic equation of state 
p = KrC and use the first law of thermodynamics for 
isentropic flows, d{p/n) -f pd(l/n) = Tds, (s is the en- 
tropy per baryon and T the temperature), we find 



7-1' 



(8) 



Now, assuming that the flow is relativistic we can ne- 
glect the rest mass energy density po with respect to the 
internal energy density Kn* / {'^ — 1), i.e., the flow satis- 
fies a barotropic equation of state, p = (7 — In this 
case, the baryon number density can be written in terms 
of the enthalpy as 



where k = 
written as 



( lui 



n = , (9) 

. In this case, equation ^ can be 



7 



7- 



1 



ln/i($)],„$'" = 0. 



(10) 



The simplest case is found when 7 = 2, Eq. ([7]) reduces 
to a linear equation. In this case the barotropic sound 



light, i.e., we have a stiff equation of state. In the general 
case (1 < 7 < 2), the differential equation is nonlinear 
and the barotropic sound speed, = (7 — 1), is less than 
the speed of light. 

The normalization condition, UaU°' = — 1, gives us a 
relation between the enthalpy and the scalar field. 



h = v/-«',a$^". 



(11) 



This relation will be also useful to determine the baryon 
number density. 

In this section we have followed the work of Moncrief 
[l5| which is in accord with Tabensky and Taub [ll[ who 
consider a constant barotropic sound speed, = (7 — 1). 



III. POTENTIAL FLOWS IN 
REISSNER-NORDSTROM-(ANTI) DE SITTER 
BLACK HOLES AND RIGID SPHERES 

A. The metric 

The Reissner-Norstr6m-(anti) de Sitter metric can be 
written as 



(12) 



where / = 1 - 2M/r + Q^/r'^ - Ar^/S, M is the mass, Q 
is the electric charge and A is the cosmological constant. 
The metric (|12p has different number of horizons given by 
the roots of / = 0, which depend on M, Q and A. If A > 
(de Sitter spacetime) we may obtain four different real 
roots. For simplicity, we have taken our parameters such 
this is always true. Three of the roots are then identified 
as the inner (or Cauchy) horizon (r^), the event (or outer) 
horizon (r/i) and the cosmological horizon (rc). These 
horizons satisfy the relation < < rc- The remaining 
root is negative and no physical meaning is attached. 
Following the classification due to Brill and Hayward [l3] , 
we are dealing with a generic black hole. If A < (anti de 
Sitter spacetime), we obtain two real roots representing 
the inner horizon and the event horizon. In this case 
there is no cosmological horizon and the function /(r) 
diverges as the radial coordinate r tends to infinity. 



B. Numerical method 



speed, defined as dp /dp 



is equal to the speed of 



To solve equation (|10p in the space-time with metric 
we assume: (a) That fluid is stationary, i.e., the 
function <I> depends on time only through the addition of 
—at, where a is a constant related to the zeroth compo- 
nent of the velocity, (b) That due to the axial symmetry 
of the metric the potential $ does not depend on the 
variable tp, and (c) That the fluid is a test fluid, i.e. the 
metric does not evolve and it is given a priori. With the 
above mentioned assumptions Eq. pO)) reduces to an el- 
liptical differential equation with an inner boundary con- 
dition near the black hole or rigid sphere and an external 
boundary condition (asymptotic condition). 
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To obtain the numerical solution of Eq. (fTO)) we use 
a computational code based on a relaxation method 
with a second-order precision, five-point, finite differ- 
ence scheme. The relaxation method employed is the 
Stone's strongly implicit procedure applied in a numeric 
grid evenly spaced in both r and 9. The whole method 
is detailed in reference jl4| . 

In the case of a rigid sphere, for the inner boundary 
condition, the usual condition of zero normal velocity in 
the surface of the sphere is employed. In general, the 
fluid velocity must be equal to the corresponding com- 
ponent of the velocity of the surface. Since usual stars 
have gaseous surfaces (not hard), this condition does not 
describe a typical flow around a star. In special astro- 
physical situations, like relativistic flows around a neu- 
tron star, this condition can be valid. For a black hole we 
use the condition that the fluid particle number density 
must remain finite on the black hole horizon 0. This 
leads to the numerical condition ^] that near r = r^ (the 
black hole horizon). 



where dr* = dr/ f defines the tortoise radial coordinate 
p^. With the assumption that = —o, < 0, the con- 
dition (|13p can also be satisfied in our case and it will 
be taken as the inner boundary condition. For numerical 
applications of this condition see d, 0] ■ 

On the outer boundary we apply the dirichlet condi- 
tion imposing that the fluid is homogeneous and has a 
constant velocity parallel to the z-axis of the coordinate 
system, say 

<P = -U^t + Uaorcos9, (14) 

where = ([/^,Uoo) = {I - vlr^^\l,v^). This 
condition is applied in a region located usually between 
the surface of the sphere (or black hole horizon) and 60 
Schwarzschild radius. Note that the presence of a cos- 
mological horizon (for A > 0) does not represent any 
further difficulty. Our domain of interest is delimited to 
the asymptotic boundary condition, which we carefully 
apply in a radius much smaller than re- 

To solve the non-linear equation (jlOp we first put the 
non-linear term equal to zero and calculate the linear 
part of the equation to find <&. This solution is used as 
an initial guess for the non-linear problem. Then, from 
pip we compute the fluid enthalpy. Finally, with this 
information we compute the non-linear term that is in- 
troduced in the non-linear equation to flnd new values of 
<f>. The process is repeated until the sum of the fractional 
change in the enthalpy for all interior points of the grid 
in one iteration is less than a certain error ^. Usually 
less than 10 iterations were required to reach an error of 
^ < 10~^. Obviously, for the linear case we need only 
one iteration. The number of iterations depends on the 
value of 7. We found that for values of 7 < 5/3 our code 
does not converge, this is because the enthalpy becomes 



negative in some points of our grid. So, this value of 7 
may represent a physical limit to our model. Moreover, 
when 7 = 1 the fluid is pressureless (dust) and the fluid 
flow is geodesic, i.e., no longer obeys pU)l . 

The code was tested for the case 7 = 2 using the ana- 
lytical results for the steady flow of a fluid in the presence 
of either a hard sphere [1] or a black hole P, . In this 
case the numerical solution agreed with the exact solu- 
tion within an error better than 1% for the radial velocity 
and angular velocity. This accuracy is sufficient for our 
purposes. 

Another aspect of the fiuid dynamics is the accretion 
rate of matter into a black hole. This accretion rate can 
be computed from ([5]) and ([5]), we find 

N = - [ nW^/^dS, ^ ~ [ Kh^^^rg'^^'./^dfl, 
Js Js 

(15) 

where the integration is performed in a two-surface 
sphere centered on the black hole. The exact form of (|15p . 
for 7 = 2, is known [l|,[2|- It is iV = Anrf^nooU^, where r/j 
is the event horizon and rioo and are the asymptotic 
density and zeroth component of the 4-velocity, respec- 
tively. The quantity rioo is a constant in this case. We 
test our code to compute the accretion rate with this 
exact solution. We found an error less than 1%. 



IV. NUMERICAL RESULTS 
A. Black holes 

In Fig.[T] we show the streamlines and density contours 
(n/n) for a potential flow with 7 = 2 and 7 = 1.7 through 
a Reissner-Nordstrom black hole for different values of 
the factor AM^. For increasing values of AM^, and flxed 
Q and 7, we see, from the upper and middle graphs of 
Fig [U that the density in the incident region increases 
and in the opposite side decreases. On the other hand, 
for decreasing values of AAf ^, we obtain, from the middle 
and lower graphs of Fig. [1] that the density decreases in 
the incident region and increases in the opposite side. 
These behaviors can be seen as a shift of the density 
contours (more evident for n/n = 2.0) and an expansion 
or contraction of the density contours (more evident for 
u/k = 1.5). These effects are enhanced for lower values 
of 7, e.g. the contour n/n = 1.5 does not even appear 
in the upper graph when 7 — 1.7. It must be noted, 
however, that the adopted values for AM^ are choose in 
order to magnify their effects on the flow. 

In Fig. [1] we use a flxed value for the electric charge, 
say Q — 0.9M. Further simulations show that the quali- 
tative properties of the density contours are independent 
whether the black hole has electric charge or not. Never- 
theless, the density contour values are more sensitive to 
changes in A if the black hole has small or null electric 
charge. In general, all the effects are more evident for 
lower values of 7. 
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FIG. 1: Numerical results for the streamlines and density contours (n/n) for a stiff fluid (7 = 2) and a fluid with 7 = 1.7 
flowing through a black hole with an electric charge Q = 0.9M. Three values for AM^ were considered for comparison (null, 
positive and negative cosmological constant). The fluid flows from left to right with Voo =0.6. 



Another quantity of interest for flows through a black 
hole is the accretion rate of matter given by Eq. ([15]). 
We numerically integrate Eq. (fT5|l . using Simpson's rule, 
to find that the presence of a positive or negative cosmo- 
logical constant leads to an increase or decrease in the 
accretion rate, respectively. Furthermore, the increase 
due to A > is numerically equivalent to the decrease 
due to A < 0. The accretion rate is slightly larger if the 
black hole has small or null electric charge. 

Regardless the value of A and Q, the accretion rate 
always increases for lower values of 7. Also the accretion 
rate is more sensitive to variations in both Q and A if 7 
is less than 2. For example, if Q = 0.9 M, the influence 
of the cosmological constant is about ten times greater 
when we lower the value of 7 from 2 to 1.7, and for Q = 



is almost twenty times greater. Moreover, examining the 
streamlines pattern, we note that, if we increase or de- 
crease the value of A we also increase or decrease the 
number of streamlines that fall into the black hole, re- 
spectively, see Fig [1] For lower values of the parameter 
7 the streamlines flux through the black hole increases, 
this behavior is in accordance with the previous numeri- 
cal results for the accretion rate. 



B. Rigid spheres 

In Fig. [21 we show the streamlines and density con- 
tours {11/ k) for a potential flow with 7 = 2 and 7 = 1.7 
through a Reissner-Nordstrom rigid sphere for different 
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FIG. 2: Numerical results for the streamlines and density contours (n/n) for a stiff fluid (7 = 2) and a fluid with 7 = 1.7 
flowing through a hard sphere of radius R = 2.5M with an electric charge Q = 0.9M. Two values are considered for AM^ to 
contrast to the case of A = 0. The fluid flows from left to right with Voo = 0.6. 



values of the factor AA'P . We see that for increasing val- 
ues of AM^, with fixed Q and 7, the density contours 
expand, which is associated to an increase in the density 
around the rigid sphere. Also, the effect of the cosmolog- 
ical constant is to change the oblateness of the density 
contours. For large values of A the contours are more 
spherical, this is easily noted from the u/k = 1.9 den- 
sity contours of the 7 = 2 column graphs. When AM^ 
decreases, we obtain the opposite behavior. These ex- 
pansions and contractions of the density contours are en- 
hanced when we lower the value of the parameter 7, this 
can be seen through a direct comparison from the graphs 
of the different columns of Fig[2l For all parameters, the 
streamlines suffer little changes. 

In addition, further numerical simulations show us 



that, for a fixed 7, the density of fluid around a hard 
sphere is more sensitive to changes in A if the rigid sphere 
has a large electric charge. This is contrary to what hap- 
pens in the black hole case where we find that the effects 
are less sensitive when it has a large electric charge. 



V. CONCLUSIONS 

We have numerically analyzed the effect of the elec- 
tric charge and both, positive and negative, cosmological 
constants on a potential flow of an ideal fluid with poly- 
tropic equation of state p — Kn^ around a black hole 
and a hard sphere by means of the streamlines, the fluid 
density contours and the accretion rate (black hole case 



6 



TABLE I: Qualitative description of the influence of the parameters {Q, A, 7) and their combined changes on the potential flow 
through black holes and rigid spheres. 

Variable Variation performed Summary of effects: Black holes Summary of effects: Rigid spheres 

parameter 



Increase of the ratio 
Q'^/M^, with A = 0. 



Decrease of A'^. Increase of the density 
around the black hole. Expansion of the 
density contours. 



Decrease of the density around the 
sphere. Contraction of the density con- 
tours. 



A 



Increase/Decrease 
AM^ and Q = 0. 



of Increase/Decrease of A^. In- 
crease/Decrease of the density over 
the region of incidence of the fluid 
and decrease/increase of the density 
on the opposite region. Slight ex- 
pansion/contraction and shift to the 
left/right of the density contours. 
Increase/Decrease of the number of 
streamlines that fall into the black hole. 



Increase/Decrease of the density around 
the sphere, mainly over the region of 
6 = 7r/2. Expansion/Contraction of the 
density contours. 



Q,A Increase/Decrease of 

AM^ and Q 7^ 0. 



7' Decrease of 7 to values 

smaller than 2. 



The effects due to the presence of A are 
qualitative the same as those observed 
if Q = but less evident. 

Increase of N. Increase of the density 
around the black hole. Expansion of the 
density contours. Increase the number 
of streamlines that fall into the black 
hole. The effects due to variations of Q 
and/or A are enhanced. 



The effects due to the presence of A are 
qualitative the same as those observed 
if Q = but more evident. 

Increase of the density around the 
sphere. Expansion of the density con- 
tours. The effects due to variations of 
Q and/or A are enhanced. 



"This case has analytical solutions for 7 = 2 

''The effects of 7 < 2 had been treated in other physical situations. Sec for example references Q] 



only). The main effects observed in our numerical simu- 
lations are summarized in Table IH 

As a final comment, we must say that results in cosmol- 
ogy indicates with 99% of confident that A > with 
possible values of A w IQ-^Vm^ [H, [H, [13]. A value 
such as AM^ « 10^'' implies in an object with M « lO""^^ 
solar masses. If we considered that the supposed super 
massive black holes in the center of some galaxies have 
masses of order M « 10~^, we conclude that the detec- 



tion of the influence of the cosmological constant near 
such objects is going to be difficult. 
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